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SUMMARY 



The relation between, the elevator hlnge-moraent param- 
eters and the control forces £oi^ changes In fox*ward speed 
and In maneuvers Is shown for several values of static 
stability and elevator mass balance. 

The stability of the short-period oscillations Is 
shown as a series of boundaries giving the Limits of the 
stable region in terms of the elevator hinge-moment param*- 
eters. The effects of static stability, elevator moment 
of Inertia, elevator mass unbalance, and airplane density 
are also considered. Dyn^ic instability is likely to 
occur if there is mass unbalance of .the elevator control 
system combined with a small restoring tendency " (high 
aerodynamic balance). This inatabillty can be prevented 
by a rearrangement of the unbalancing weights which, how- 
ever^ Involves .an Inci^ease of the amount of weight neces- 
sary. It can also be prevented by the addition of 
viscqu? friction to. the elevator control system provided 
the .alrpland center of gra^vlty is not behind a certain 
critical position. 

Poip high values of the density parameter^ which oorre-. 
spond to high'altlt\2des 'of flight,, the addition of moderate' 
amounts of viscous friction niay be destabilizing " even when 
the airplane is -statically stable. * ■ tri this case, .in- 
creasing the viscous friction makes the oscillation stable 
agaln« The condition In which viscous" friction causes 
dynamic, injaitablllty of a statically stable airplane is 
limited to a definite ^angis of hlhge -moment parameters. 
It Is shown that, when viscous friction causes increasing 
oscillations, solid friction will produce steady 0801118-* 
tlons having an amplitude prdportlonal to .the amount of 
friction. 
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INTRODIJCTIOII ■ 



-The effects of aerodTnamic balance and mass imbalance 
of the elevator on the dynamic stability of the airplane 
are discussed In a previous report on control-free sta- 
bility (reference l") • It was found theoretically In ref- 
erence 1 and verified In flight (reference 2) that, If the ■ 
elevator Is too closely balanced aerodynamlcally and has 
a sufficient amount of mass unbalance (which tends to de- 
press the elevator). Increasing oscillations of short 
period may occur. Other flight tests (reference 3) 
showed, however, that mass unbalance of the elevator con- 
trol system Improves the static stability of an airplane, 
that Is, Increases the slope of the curve of stick force 
against speed- In level flight and of the curve of stick force 
against normal acceleration In maneuvers. .Subsequent 
work (reference has Indicated that a control siirface 
with positive floating tendency (tendency to float against 
the relative wind), when used as a rudder. Is effective 
In Improving control-free static stability. A theoreti- 
cal analysis (reference 5) showed that a rudder having a 
positive floating ratio may, under the Influence of s.olld 
friction In the control system, build up steady oscilla- 
tions of the airplane and rudder. These steady oscilla- 
tions have been observed in flighty tests (reference 6). 
These results suggested an investigation of the behavior 
of an airplane equipped with an elevator having a positive 
floating tendency. This type of elevator was not con- 
sidered In any of the previous investigations. 

The purpose of the present report Is to make a 
theoretical analysis of the control-free- longitudinal 
stability of an airplane, which takes account of" this 
current trend toward a positive floating tendency in 
control-surface design and covers. In general, a much wider 
range of parameters than the investigation of reference 1> 
These parameters Include, for the elevator control system^ 
restoring tendency, floating tendency, mass unbalance 
(bobwelght control), moment of inertia, and viscous and 
solid friction .and, for the airplane, density and center-of- 
gravity position. 

The method of analysis of dynamic stability is based 
on the classical theory of .Bryan and Balrstow extended to 
Include movements of the controls and their couplings 
with the airplane motions. Friction is treated In the 
same way as in the approximate method of reference 3f 
which solid friction Is replaced by an equivalent viscous 
friction. 



Before the analysis of dynamic st^bLllty Is presented, 
some dlsousslon Is given of the effect of the various 
parameters on the elevator forces for trl^ and-' for aocel«* 
eratlon ^ characteristics considered Important to flying 
qualities, The stability of the- short-period oscilla- 
tions,* with and without friction In the control system. 
Is then considered* The effects of wel^ts added to the 
system to modify the static and dynamic stability are dis- 
cussed. The trends to be- ejcpected are Illustrated by a 
series of calculations and charts based on a typical air* 



plane. The stability of the long-period (phugold) 
oscillations Is not discussed because of Its relative im« 
Importance • 

SYfJBOLS 

Aw wing aspect ratio 

At tall aspect ratio 

A,B,C,E,P coefficients In stability equation 

b wing span 

Ch elevator hinge-moment coefficient 3 ^\ 

Ch* frlctlonal hinge-moment coefficient/^ ^ — ^ • ^ 

Cho applied hinge-monent coefficient 

Cl airplane lift coefficient (M^) 

lift coefficient of tall 

Cm pltchlng-moment coefficient about airplane 

center of gravity 

c wing chord 

Ce elevator chord 

D differential operator (d/ds) 

' constant tern In vtaMllty equaCtlon 



stick force 

stick-force gradient In maneuvers 

stick- force gradient for level fl; 

acceleration of gravity 

hinge moment; positive when tends to depress trailing 
edge 

mass moment of elevator about Its hinge; positive 
when tallheavy 

mass moment of control stick about Its pivot; 
positive when stick tends to move forward 

frlctlonal hinge moment 

He + rHa 

rhg + he 

4He/pSeC©c 

ij-Hg/pSeCeC 

moment of inertia of elevator about its hinge 

moment of inertia of control stick -about Its pivot 

le - rls 

le + r^ig 

8le/pSeCec2 

8ls/pSeCec2 

radius of gyration of airplane about Y-axis 

2ky/c 

distance between airplane center of gravity and 
elevator hinge 

2Lh/c 

length of control stick 




H pltohlng moinont about airplane center of gravity 
p mas a of airplane 

o Ni/2 number of cycles required for oaclllatlon to daiop 

^ to half unplltude 

n normal acoeleratlon par g of alrplaxie dufe to 

ourvatxire of flight path; aocelerometer reading 
minus component of gravity force ■ 

P period of oaclllatlon, seconds 

q dynamlo pressure 

Sq elevator area' 

St tall area 

wing area 

s distance In half-chords (2Vt/c) 

/p time required for oscillation to damp to half 
amplitude, seconds 

t time 

u = Av/y 

V forward velocity 

AV change In foiTward velocity from trimmed value 
W weight of airplane 

X longitudinal force; positl'Ve forward 

Xa.c. distance of center of gravity from aerodynamic 

center; positive when center of gravity is ahead 
of aerodynamic center 

Z normal force; positive downward 

a angle of attack 

at angle of attack at tall 
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deflection of elevator; positive for downward 
mptlon of trailing edge 

7? amplitude of elevator oscillation 

€ angle of downwash 

r control gearing (Q^/Sq) 

e angle of pitch of airplane 

&g deflection of control stick; positive for forward 
motion of stick 

\ complex root of stability equation 

g,ri real and Imaginary parts, respectively, of \ 

\jL airplane -density parameter (m/pS^b) 

p mass density of air 

Whenever |x, V, a, at, 9, 6, Da, D9, D6, 

and D^a are used as subscripts, a derivative Is indi- 
cated. For example, Xy = dx/6V and Chj^g = ^0^/^05^ 

Whenever a dot Is used above a symbol. It denotes dif- 
ferentiation with respect to time. 

All angles are measured In radians. 



METHOD OF ANALYSIS 



Four degrees of freedom - forward speed, angle of 
attack, angle of pitch, and elevator deflection - are 
generally Involved In the problem of control-free sta- 
bility. To each degree of freedom, there corresponds an 
equation of equilibrium between Inertlal and aerodynamic 
forces or momenta. By use of wind axes, the four equa- 
tions become, for level flight. 



XV + Xaa + ZqG " n 

AV + ZjjO - aV(d - 8) 

E7 A7 + Eaa + B^a * Byfi * Hgd + E66e + Hsfie " + + He(Va - T)6 -f 

♦ rjlg(r(5e - ^ + HB(Va - 



which oan be vritten in nondimensicmal form as 

Cl 
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(Xy + 2A^)\x + Xgo. + -5-0 - 0 

- 0 



°V * * ' * Va°^" * ' * (^^6 * V ' ^2^- " ° . 

In applying equations (1) to dynamic stability, certain approximations may be made. For 
instance, 8ho3rt-period oscillations (of the order of 1 sec) involve negligible changes in forward 
speed, ivhich may thereforo be neglected in studying the short-period oscillations. In foot, the 
period and damping of these oscillations can be obtained to a high degree of accuracy by Using 
only the last three of the equations (1) and setting u " 0. 



00 



Equations (1) then beoone 

(^a * V * "^T^a^y + (Cmj^ - 2A^2D)De * (Cms * C^j^X " 0> (2) 

fch^ + (ChDo " *')° * ^Vo^** * P^D9 * ^ "(*^^^ * il)D[De + (Chg + ChjjgD - i2D2)6e - 0 



»•.■ By setting 

a = Ooe^^t D6 - (D8)o«?^t - 6j»\t 

It oan be shown (refsrenoe 7) "U^t \ must be a root of a quartic equation fomed by writing 

— + 2A^\ -2A^ 0 

The resulting stability equation may be written as 

+ B\3 + CX? + E\ + P - 0 
where A, B, C, B, and F are funotions of the stability derivatl-res. 



(3) 



Hhe study of the effects of different parameters on 
the control-free stability was made by a series of compu- 
tations for an average airplane having the characteristics' 
given hereinafter. T!h.e current trend toward a positive 
floating tendency In control-surface design suggests the 
USB of Ch^^ and Chg as the fundamental variables to be 

used In expressing stability and control characteristics. 
The results are presented as a series of figures that show 
the relations between Ch^^ and ChQ which, with the 

t 

other derivatives flxed^ satisfy the conditions fbr neutral 
dynamic stability and neutral static stability. 

A curve for neutral dynamic stability Is the boundary 
dividing the region of Increasing oscillations from the 
region of damped oscillations and Is obtained 'from Routh*a 
discriminant 

BCE - AE2 - Pb2-=0 
The condition for neutral static stability Is that 

P = 0 



The stability equation (3) has four roots. A pair 
of complex roots Indicates an oscillatory mode and a real 
root Indicates an aperiodic mode. The real part of the 
complex z^oot determines the damping; the Imaginary part 
dete]?mlnes the period of the oscillations. More specifi- 
cally. If there Is a pair of complex roots ' 

^ = S. ± ill 
the period In seconds Is given by 

^ " 2V T] 

and the time In seconds to damp to half amplitude Is given 

^1/2 - - 2V £^ 

For an alrplGine at constant speed, there may be two 
oscillatory modes, there may be only one oscillatory mode, 
or the motion may be entirely aperiodic. In cases In 
which there are oscillatory components, one of the oscilla- 
tions may be poorly damped and even become unstable. 
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The average airplane on wMch the calciilatlons of 
this report are based Is of conventional design. !l*he 
characteristics of the airplane are 



% 1.5 

Lh/o = lh/2 3.3 

r 1 

St/Sw 0.18 

la' " 2 

Se/St 0.55 

At 1+.5 

The. basic stability derivatives and parameters obtained 
from these airplane characteristics by methods shown In 
appendix A are 

% ^-5 c^De -^5.3 

°Ltat Crnoa ::::::: • -8-9 

O.W Chi3^ 3-22Chat 

Cme -i-5ij- Cmj^a^t 25.2 

°mD6 -0.97 Ch^Sa • • • -lO^SSCha^ 

Chjjg (with no friction). -1 

The following parameters of the airplane were varied? 

\i airplane -density parameter 

Cjn^ control-fixed static-stability parameter 

The following parameters of the elevator control system 
were varied: 

Gil floating tendency 
"-t 

restoring tendency 
CjiDs e leva tor- damping parameter 
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le moment-of-lnertla parameter of elevator about Its 
hinge 

I3 moment-of- inertia parameter of control stick about 
Its hinge 

h mass-moment parameter of elevator control system 
about elevator hinge 

h^ mass^moment parameter of elevator alone about its 
hinge 

As has been pointed out^ the stability boundaries 
were plotted, in most cases, in terms of ^h^^ Ghg 

as the variables of the coordinate system. In analyzing 
the effects of friction in the control system, ChQ 
and Gh]35 ^^^^ used as the plotted variables in some fig- 
ures whereas Cha^ ^hs were used in others. The 

effect of the other parameters is found by varying theiji 
one at a time, through a range of values, and showing for 
each parameter a series of stability boimdaries. 

The size of the airplane, together wlth'wlng loading 
and altitude, are combined in the parameter \i, which 

is J^- . A variation in u. thus could be due to a 
pSwb 

variation in size, wing loading, or altitude, or any com- 
bination of these. The range of valioes of [i covered 
in the present report together with some typical corre*"- 
spending values of wing loading, altitude, and size are 
given in the following table; 





wing loading 
(Ib/aq ft) 


Altitude 
(ft) 


Mean wing chord 
(ft) 


ij-.l? 
12.5 
37.5 


ho 

^0 

h.0 


Sea level 
Sea level 
35,000 


21 
7 
7 



The range of CniQ and the corresponding center-of -gravity 
positions are as follows: 





. c . • 

(fraction of mean wing chord) 


-0.352 


-0.05 


0 


0 


.252 


.05 
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The ranges of values of the other parameters, for a small 
airplane (chord, 7 ft), are as follows: 



1 ■. 


Moment of Inertia of elevator control system 

(alug-ft2) 


0 


0 


s 


1.6 


ii. 


5.2 



Stick force 
(lb) 



a 


At sea level 


At 35,000 ft 


0 


0 


0 


10 


57 







dH/d6 




V 


(lb-ft/100 raph/deg/aoo) 




At sea level 


At 55,000 ft 


-10 




1.95 


-100 


58.5 


19.5 



STATIC STABILITY AND RELATION TO CONTROL FORCES 



The connection between the static stability and the 
airplane and control parameters Is established to assist 
In the Interpretation of the resiolta obtained herein- 
after. Equations (1) can be applied to static stability 
by setting all terms containing D and 1)2 equal to 
zero and solving the resulting equations simultaneously 
for the variation In foi'ward speed with an applied 
elevator hinge moment. For level flight, 6 Is filso 
zero and the resulting equations are 



Solving give a 



variation of stick force with fractional change In 
forward speed Is 



Is^L 



If eTfects of allpatream on the tall are neglected, 

KCt 

Cm =0. A3 shown In anpendlx A, C^^ " — . In- 

aertlng these values in the expression for Cho/^ shows 
that Fu Is Independent of forward speed. 

IJtxQ variation of control force with normal acceler- 
ation In a steady pull-up, with no change assumed In for- 
ward speed (see reference 8), can be found from equa- 
tions (2) by equating to aero all terms containing D 
except D9. This procedure implies that the normal ac- 
celeration is due entirely to curvatuire of the flight 
path LBe The equations become, for an applied hinge 
moment, 

- ZA^ii De =0 

Cma°- ■*• CnoeD-Q * CmgOe = 0 



frcaa which 

Ohp _ -^*i#ChaCinfl- ^Lg^' ^e^Lg^S* M^ChsSi* 
If the normal acceleration is ng, 

1 ^ 



Pn = 



Cho pSeCeCg 
Be IJT^r 

These fomulas for P^^ and Pn are equivalent to equa- 
tion (1) of reference 9 and equations (27) and (28) of 
reference 10. 

The formulas- indicate that the stick- force gradi- 
ents and Pj^ are dependent on most of the afore- 
mentioned airplane and elevator parameters. Figures 1 
to 5 show the variation of these stick- force gradients 
with the parameters ^h^^* ^^q' ^"^a' ^^^^ 

The gradients are Independent of speed, although only 
within the limits of the asaumptlcna rade in the analysis, 
namely, neglect of power and of compresslb.lllty effects. 
The gradient P^ can be used to get the stick force i*or 

only a small change in forward speed because the stick 
force is not directly proportional to the change in speed. 
The stick force in a steady pull-up Pj^, however, is 
proportional to the normal acceleration provided the 
control deflection is not so groat that "the basic assumTD- 
tlon of linearity is violated. 



The line P-u " 0 is the boundary for true static 
stability - that is, = 0 Is the condition for zero 
variation In stick force with forward speed in steady 
flight. This condition is the same as that obtained by 
setting P = 0, where P la the constant term of the 
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bh-order stability equation obtained from equations (1). 
subsequent figures it is ^called the divergence boundary. 



slxth-c 
On SI 

ffltie line Pn = 0 is the boundary for apparent static 
(or maneuvering) stability and is the condition for zero 
variation in stick force in a steady pull-up. On the 
unstable side of " 0» ^ slow divergence occxars that 
Is noticed by the pilot as an unstable variation o;f 
stick force with forward speed. ?3ie stick force due to 
normal acceleration in a ptill-i5) Is stable, however, 
unless the conditions are such that the airplane is 
operating on the unstable side of Fn ^ 0- 

Piguros 1 to 5 Indicate that the parameters have 
the same effect on and Pn except that the altitude 

affects only P^. They show that the stick-force gradi- 
ents .on an airplane of given tall size and center-of- 
gravlty position may be increased by making the floating 
tendency Cha^ more positive or by mass unbalancing 

the elevator control system to depress tho elevator 
(make it tallhoavy). The effect of the restoring tend- 
ency Ciig on the stick-force gradients depends on the 

relative position of the center of gravity and the 
aerodynamic center. If the center of gravity is ahead 
of the aerodynamic center (airplane stable with controls 
fixed), increasing the magnitude of Chg increases the 

stick-forcG gradients. If the center of gravity Is 
behind the aerodynamic center, this effect on Fyx 
reversed; the effect on is not reversed, however, 

until the center of gravity is well behind tho aero- 
dynamic center (in this case, about O.O^c at sea level 
axui 0.020 at 30,000 ft). If Ou ^ 0, the stick forces 



are Independent of the position of the airplane center 
of gravity. 

Increase in altitude will either Increase or decrease 
Pn, depending on the hinge-moment parometers, The solid 
line in figure 5 is the locus of values, of Ch^^ and Oh^ 

for which Fj^ is Independent of altitude. For points 
to the left of this linsi decreases with altitude; 

for points to the right of this line, P^^ increases with 
altitude. Olhls line is determined by the relation 
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wMch, for the case of figure ^, becomes 

Cha^ = 1.50 ChQ 

Another method of Increasing the stick-force gradi- 
ent In level flight consists in applying a constant 
hinge moment to the elevator by means of a spring or 
bungee. The effect of the spring on the gradient F\i 
is due to the derivative Gh^ which depends in the same 
way on the constant hinge moment, whether it is caused by 
a weight or by a spring, A bungee^ which tends to de- 
press the elevator, will therefore Increase the stick- 
force gradient in level flight P-u» The effect of the 
bungee on the stick-force gradient in accelerated flight 
Pn will be zero because its action depends solely oii 
• changes in forward speed. Its eff&ct on the short- 
period oscillations will be zero for the same reason. ■ 



DYWAMIG STABILITY 
ITo Friction in Control System 

The stability of tlie short-period oscillations with- 
out frlct?on is shown in figures 6 to 11, which also show 
the boundaries for true static stability (divergence 
boundaries). Pigure 6 is an example of a more nearly 
complete presentation of the stability data than subse- 
quent figures because it shows the variation of damping 
and period of oscillation with the hinge-moment param- 
eters Ghfl and Chg for certain fixed values of tho 
t 

other parameters. The damping, which is proportional 

to £, increases with the magnitude of Ch5* ^e period, 

proportional to ^, decreases as GhQ^ increases. 

Another way of presenting this additional stability data 
is shown in figure 7* which gives the number of cycles 
the oscillation performs before it damps to half fiimnll- 
tude. It is clear from figure 7 that the oscillation 
is very well damped unless tho restoring tendency is 
close to zero. In this particular case, only one oscil- 
latory mode exists. Inasmuch as there are only three 
roots in this case (because I2 and ii = 0), the other 
root is always real and Is of no particular significance 
in dynamic stability. In cases in which an additional 
oscillatory mode exists. It is. always stable. 
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. The effeot of oenter-of -gravity position on the 
atahlllty bf the short-period osclllatl'ona Is shown In 
figure 8. The shift In the dynaialc-stablllty boundary, 
for the comparatively large shift In center of gravity 
shown. Is practically negligible. Many of the subse* 
quent figures, in which zero static stability Is assumed 
to facilitate computation, therefore are valid for air- 
planes having a margin of static stability. 

The effeot of moment of inertia of the elevator 
control system on the dynamic stability is shown in fig- 
ure 9, which gives typical values of the moment of 
inertia. The effect is slightly destabilizing especially 
for high values of . The destabilizing effect of 

the moraont of Inertia of the elevator Is greater than that 
of the moment of inertia of the control stick. Because 
the accuracy gained by including moment of inertia is 
small compared with the saving In labor gained by neg- 
lecting It, moment of Inertia of the elevator control 
system was set equal to zero in the subsequent calcula- 
tion. As a result, the stability equation becomes a 
cubic and subsequent figures are somewhat unconserva- 
tlve. 

The effect of density parameter \i on the dynamic 
stability is shown in figure 10. Increase of \i has a 
slight destabilizing effect. 

As has been shown, mass unbalance of the elevator 
control system improves the static stability (fig* Ij.) . 
The effect on dynamic stability Is unfavorable, however, 
as shown in figure 11, The value of mass unbalance 
shown corresponds to a bobwelght that is located at the 
airplane center of gravity and requires a balancing pull 
of 37 pounds on the control stick of a pursuit airplane 
at sea level. Increasing oscillations occur if the 
aerodynamic balance is too high (low magnitudes of Chs)^ 
especially for negative values of C-^m 

The effect of the location of the bobwelght is shown 
in figure 12. Each curve represents a different arrange- 
ment of bobweli^ts and all arrangements give the same 
stick force. The solid line corresponds to a weight at 
the airplane center of gravity (for practical purposes, 
at the control 3tick). The short-dash line corresponds 
to a weight at the elevator. - The long-dash line corre- 
sponds to two weights - one at the elevator, which tends 
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to make It noseheavy; the other at the control stick, which 
gives the elevator a sufficiently pow.erful tallheavy moment 
that the resultant stick force Is the same as with the 
single weight. In the particular case represented, the 
noseheavy moment due to. the weight at the elevator Is 
equal to the tallheavy moment due to both weights. Moving 
the single weight from the control stick to the elevator 
has a large destabilizing effect. Overbalancing the 
elevator while the stick force Is kept constant has a con- 
siderable stabilizing effect. This method of preventing 
Instability has the disadvantage, however, of Increasing 
the total amount of unbalancing weight required. . In the 
case shown, the weight Is Increased to three times Its 
original size. Another disadvantage Is the rearward shift 
In center of gravity of the whole airplane due to addi- 
tional weight at the elevator, (See airplane parameters 
given In "Method of Analysis.") The destabilizing effect 
of the increased moment of inertia associated with the 
added weights' was found to be very small, especially for 
negative floating tendency. 



Effect of Viscous Friction in Control System 

In the preceding section, a constant value of the 
elevator -damping parameter Chjjg was assumed. This 

value was due only to aerodynamic damping. The effects 
of viscous friction in the elevutor control system are 
obtained by considering Chj^j^ as an additional variable. 

This variable can be Introduced, es in the preceding sec- 
tion, by showing a series of boundaries in the ChQ-^Ciig 

plane for various values of ^hj)5« ^® general nature 

of the effect of friction is shown first, however, by 
presenting boundaries in the ChgChj^g plane with ^h^^ 

constant and some other parameter varied. This method 
of presenting stability boundaries makes it easier to 
show the effects of other parameters such as airplane 
center-of -gravity position and density when friction is 
Introduced • 

The effect of viscous friction on the dynamic sta- 
bility, for various conditions, is shown in figures 15 
and llj. for \i = 12.5 and \i = 57*5^ respectively. 
Figures 13(a) and lij-Ca) refer to the . mass-balanced 
elevator control system; figures 15(b) and lii-(b) refer to 
the tallheavy elevator control system considered in the 
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preceding section. It Is shown that. If the airplane 
oenter of gravity Is ahead of a certain point, ^ the 
Instability caused by the tmbalanced elevator can be 
removed by adding viscous friction to the control system, 
o TblB criticELL center- of -gravity position is behind the 
!B aerodynamic oenter, and its distance from the aerodynamic 
»^ center decreases as the density parameter (jl increases ^ 
Hthen the center of gravity is behind this critical posi- 
tion, vlscotia friction has a destaMllzlng effect. Xbese 
opposite effects of viscous friction are shown in the 
(^a^Ohfi plane in figures I3 and l6. lllhen the center 

of gravity is slightly ahead of this critical position, 
the effect of viscous friction depends on its magnitude 
and also on the value of C^. Tiie addition of a small 

amount of viscous friction la destabilizing but larger 
amounts are stabilizing. If the aerodynamic balance is 
sufficiently high iChQ ts 0) an! the viscous friction 

lies In a certain range, increasing oscillations will 
occur, in figure llj.{b), for example, if Xa.o, = -0.01c 
and Chg = -0.0^, the oscillations will be unstable when 

the elevator-damping parameter Is in the range from -2.5 
to -76. If Chfl is more negative than -0.086, no 

amount of elevator damping can cause increasing oscil- 
lations. As the center of gravity moves foDrward, the 
destabilizing effect of elevator daiqplng becomes less and- 
finally disappears. 

Utie effect of the densltv parameter [l can be seen 
by ccnqparing figures 13 and llj.. OSie critical center-of- 
gravity position approaches the aerodynamic oenter as \l 
increases, DiOien ^ = 12*3, elevator dainping always has 
a stabilizing effect provided x^.o. positive. When 

V' = 5$ elevator damping may be destabillaing over a 
small range of Chna QhA oven when Xo a is 

positive To. 050). ^ 

When the center of gravity is slightly ahead of the 
af ore-zoentioned critical position (which is beblnd the aero- 
dynamic center), the conditions under which elevator 
daiqping may cause dynamic Instability may be advantageously 
represented in the Ohoj-Qha - plane. If a series of sta- 
bility boundaries are drawn in that plane for various 
values of elevator damping, they will all be confined to a 

^Sinoe this report was written, this point has been 

found to be where ^ ^ 0 sometlmss oallod the stlok-flxed 
maneuver point. 
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reglQn bounded by a line that will be called the boundary 
of complete damping. An Illustration of two methods of 
constructing this boundary la given In figure 17- If -a 
series of boundaries In the Ch^ Chs plane are drawn for 

various values of the damping, the common tangent of all 
these curves Is the boundary for complete damping. ^Is 
botindary can also be drawn by plotting the minimum values 
of. C^g obtained from plots of the type shown in flg^ ; 

ures 15 and ll|. against corresponding values of Ch^^- 

^t 

The region in the ^h(j^^hg plane between the boundaries 

for complete damping and increasing oscillations Is the 
region where the addition of viscous- friction to the 
elevator control sTrstem may cause dynamic instability. 

That a boundary for complete damping cannot be shown 
for p. = 12,5 If the airplane is statically neutral or 
stable {x^^rj, is zero or positive) may be seen from 

figure 15. It is possible, however, to show a bovmdary 
for complete damping \mder these conditions for \i = 37 •5« 
Figure 18 shows these boundaries for ^Cq.c. ^ 0 and for 
the critical value x^.c. = -0,017c, for both a mass- 
balanced .elevator and a mass-unbalanced elevator. The 
boundaries for increasing oscillations and divergence are 
also shown. For the case of the mass-balanced elevator 
(h = 0), the boundary for complete damping is a straiglit 
line through the origin and therefore corresponds to a 
fixed ratio of the floating and restoring tendencies, or 
floating ratio. ^levator mass unbalance decreases the 
region of complete damping. 



■ Effect of Solid Friction in Elevator Control System 

.The boundary for complete damping has an Important 
bearing -on the effect of solid friction on dynamic sta- 
bility. In order to calculate ttiis effect, the solid 
friction is replaced by an equivalent viscous friction that 
would dissipate energy at the same rate. This condition 
gives an equivalent 

%o " i ^ 

for a sinusoidal motion of the elevator wltb amplitude^ & 
and angular frequency tj. 
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It can be shown that If viscous friction Is desta- 
bilizing-, as In figures l6 to iQ, solid friction may lead 
to steady oscillations having an fiunplltude proportional 
to the amotint of friction. Suppose an oscillation Is 
started with amplitude and frequency which result In an 
equivalent Chps that would cause Increasing oscilla- 
tions • Let this condition be represented by point 2 In 
figure 19, The amplitude of the oscillations would then 
Increase until the equivalent Chj^g decreased to the 

value that would result In neutrally damped oscillations, 
represented by point 3 In figure 19» If the Initial 
amplitude Is so low. that the equivalent viscous friction 
is In the stable region, as shown by point 1, the oscil- 
lations will die out completely. If the Initial ampli- 
tude Is so high that the oscillations are stable, repre- 
sented by point Ij., the amplitude will decrease until It 
reaches a constant value, when the equivalent Chjjg la 
again represented by point 5. Ttie value of 0^^^ at 

point 5 then determines the amplitude of the_steady 
oscillations. By solving formula (1;.) for 5, the 
amplitude of the steady oscillation Is obtained as 

5 = ^ — 



where rj and Ghpe are the values at point 5. This 
formula shov/s that the ai^niltude Is proportional to the 
amount of solid friction. 



The foregoing analysis shows that the region In the 
Ch^^Cj^g plane between the boundary for Increasing oscil- 
lations and the boundary for complete damping Is the 
region where steady oscillations may occiju? because of the 
Influence of solid friction In the control system. All 
the remarks In the preceding section concerning the 
boundary for complete damping with viscous friction conse- 
quently apply to the boundary for steady oscillations with 
solid friction. Inasmuch as these two boundaries are the 
same, within the limits of the assumptions Uivolved In the 
use of the concept of equivalent viscous friction. Steady 
oscillations due to solid friction will not occur on a 
statically neutral or stable airplane, for Instance, un- 
less IX is. very large (corresponds to a high altitude). 
Even In that case, the possibility of steady oscillations 
exists only for a comparatively small range of floating 
ratios. If the airplane Is statically unstable by a 
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sufficient amount, however, steady oscillation may exist 
over the entire range of floating ratio. 

Some calculations of the amplitude of the steady 
oscillations, expressed in terms of normal acceleration 
per unit frictional force as felt at the control stick, 
were made by the method of appendix B. The results are 
presented in figure 20, which shows lines of constant 
amplitude In the Cha^^hg pl^ne for an airplane with 

the center of gravity at the critical position referred 
to in the preceding section. Steady oscillations will 
therefore occur throughout the entire region where 
stability exists in the absence of friction. The ampli- 
tude is smallest for high values of Oix^^^ combined with 

high values of Chg* 
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The stick-free static stability of a conventional 
airplane may be Improved by making the elevator floating 
tendency more positive or by mass-unbalancing the elevator 
control system to make the elevator tailheavy. In- 
creasing the restoring tendency also has a favorable ef- 
fect provided the airplane is stable with stick fixed. 
If the restoring tendency is zero, the stick-free static 
stability is independent of airplane center-of-gravity 
position* 

The dynamic stability v/lth frlctionless controls 
depends chiefly on the restoring tendency Chg and on 

the mass balance of the elevator control system. If tlie 
elevator control system is mass unbalanced (elevator made 
tailheavy) by an offset weight at the control stick and 
if the restoring tendency is too low, increasing short- 
period oscillations nay result. This condition can be 
remedied by the use of two additional weights - one at the 
elevator making it noseheavy, the other at the control 
stick making the elevator sufficiently tailheavy that the 
combined effect gives the elevator the desired amount of 
tallheaviness. 

The addition of viscous friction to the control 
system will prevent dynamic instability provided the air- 
plane center of gravity is forward of a critical position, 
which is behind the aerodynamic center and approaches it 
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as the value of the density parameter \i Increases. If 
the airplane center of gravity Is behind tft^ critical 
position, viscous friction will have a desteihillzlng ef« 
feet, no matter -what the hinge-moment paramoters are. 
If the center of gravity is slightly ahead Qf the criti- 
cal position, viscous friction may he destabilizing for a 
limited- range -of valuejs of. viscous friction and the hinges- 
moment parameters. A low restoring tendency and a high 
positive floating tendency will ten^ .:to cauee this dy- 
namic instability. When \i Is very largb (hijjh alti- 
tude), this condition of steady osclllationc can occur 
even if the center- of .gravity is ahead of the aerodynamic 
center. 

The presence of solid friction in the control system 
may cause ahort-period steady- asclllations under the 
conditions for which viscous friction is deatabllizing. 
The amplitude of the osp-illatlons is proporibional to the 
amount of friction present. " *" 



Langley Memorial Aeronautical Laboratory, 

National. Advisory, Cpinmittee for Aeronautics, 
Langley Pield'^ ' Va.'" 
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APPEMDIX A 

■ ■• ••' .;■ r-' ■ : - . 

• ■ 36 VALUATION OP- STABIEITy-. DERIVATIVES 
• ■ - :. . . 

. •-■ •• • • • • • . . * • • 

1)9771 vatlve ^hu* ' ^^^"^^ hinge, moment - acting 
on the elevator may be expressed as 

= Chg^e.)|pv2SeOe + Hog 

• ■ ' . ' • • • •. • 

At "trim. H =.-0, therefore, •. 

••■ ■ H„g / . . 



ipV^SeCe 



■ ..Hp ,^ 2Hog . 

dH/|pV^SeCe 2Hv k^S 
6V/V ^ pVSeCe " " pV^SeCe 

Using ^pV^ClSw = mg gives 

hcgpCLSw hCL 
0. =- 

2ing 



Derivative Cma'" The parameter Cma ^7 ob- 

talned from wlnd-tuimel measurements or. If the position 
of the aerodynamic center of the complete airplane is 
known, may be calculated by the formula 
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/ o*. SA 



where 



~ = 1 - ca 

, t 

Derivatives Cigpq ^"ip^a *'^ - T^^' .v.di3a?l vaWvea 

^^™Da ^^d.: ^m])2(j arise because of -the lag between the 

•&hange in: angle of attack .at the' wing. And the'.-corre- 
a ponding do wnT/a ah at the tall. It .la a^aiined that the 
downwaah at any Instant t depends on the angle of 

• : Lh. ■''»-■..■• ; . . ■ 

attack at the Instant t '-" — \ the difference bel'ng the 

time required for the air to move from tHe v/frig t& the 
tall. If a = f (t),.. . tiilfl relation may be expressed 
as 

€ = ca - At) 

where J-^, 

it = Lh/V 

"Now ^ ' ■ » . • - " f 

■'f(t - at)" = f(t) - Af f'(tj WiiiL f"(t> - ..; *■ 

21 . 

Hence , 

« = «a 



^ - At a + -^^1^ a - . . 



■ • 2V " kV^ -3 •" 

or, because a = — Da and o = =t-s- D^a, 

C C2 



= ^a\OL - lh Da + -=—0*^0 - ... J 



and, because at = a - e , 
at = a .€ a 



la - lh Da + "^D^^ - . J 
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and 



Q.M - 0(1 - €<i) + «olh Do - «ar5- + ••• 



The part of the pltohln^-moment ooefflclent contrlb 
uted by the tall la " " 



Ih St 



at2 -Sw 

St 

at 



a(l - £0) + €dlh Da - 



The lag effectively Introduces . derivatives C^jj^* . 
Cmj^Q, . • !33ie first two of these derivatives are 



and 



= C 



Ih^ St 



Derivatives Ch^, Chpg* ^"^^ ^h-j-^^ .-- The. deriva- 
tives Ph^-f ChjjQ» aiid. Chjj2a ^-^ o^t^^'^®^ from 



Ch = Ch^i^tft 



which gives 



• Cha 

°hDa = Chat^cc^-h 
°hD2a 



27 



Derivative Cmp ^*^ ^® pitching moment jiue to 

-pitching is made up of p^rts due to pl^opeller, -^Ing, 
fuselage J and horizontal tall. The contribution of the 
tall Is by far th,e largeat and can easily be calculated. 

I ■ ■ . 

f If the airplane is rotating with. angular velocity 9, 

the increase In angle of at back at the- tall is ^hy ■» 

which results in an Increased lift on the tail given (In 
. c.p.efflcient fomj by ^ 

The resultant pitching-moment coefficient Is 

*'c -Cat V c Sw . 

and e:cT)ressing 6 as — ^D9 and — ^ as ly, gives 

c c- ■ " 

Ih^ St 
Cm = "°Lt^^-T- s; 

The conbrlbtitlon of the wing' depends on the assumed 
axis of rotation (center of gravity) but a fair average 
value will be obtained by assuming .that the center of 
gravity is at the wing quarter-chord point. This as- 
sumption gives a value 

= DB ■ ■ . 

-The total pltching-moment coefficient due to pitching 
therefore is 

: Derivative Cmfi '" derivative Cj^i^ maybe 
measured directly in a wind tunnel or may be ccmiputed 
from wind-ttuiinel data on the value of ^Lt^ 

horizontal tail by means of the form\ila 
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Deri va tl ve Clans'. " ' ^^-^ ^ ^""DO ^® 
puted from 



Ih St ct 
2 Sir c 



where ( ^^"N andl ) may be obtained from figure 1 

of reference 11, which la based on thin-wing potential -flow 
theory. 

Derivative GfiT^g -- The derivative Chpg la given 

by 

In the absence of viscous friction in the elevator con- 
trol system, the value of Cjidq may be computed from 

of reference 11 • 



where i — ) and ( -r-— 7 ^ay be obtained from figure 1 
\dD6 ^ \6D5 /3 



If a dashpot, which has a damping constant of 
K pounds per foot per second and moves a distance of Q feet 
per radian of elevator deflection, is inserted in the con- 
trol system, 

°hD5 - pVSeCeC ^^^^ 

The total value of G^j^g is the sum of equations (Al) 
and (A2). 

Derivatives Ch g^ ^^nd Che *" 'Th® derivatives Cha^ 
and Gh6 can be calculated by thin-wing*- section theory 
but the results are of doubtful aoc.uracy because of three- 
dimensional and boundary- layer effects. It is therefore 
best to obtain these derivatives from. wind-tunnel tests. 



APPEHDIX B 

■ ■ ■ 

. GALCIILATION OR NOIOL ACCSLCRATIOK'-DUS TO CSCILLATIirG ELEVATOR 

'" • . ' . • ■ . • ' • ' ■ 

The normal accoleratlon of tlio airplane^ wMch Is oqual to D(a - 0) In non- 
dimensional units^ odn-be oalculatod from 

sk-^ ' . . ■ ■ ■ ' . 

-i—. - — t. =, — 1 ^ ^ (Bl) 

8 



where xi Is the angular frequenoy of the elevator. 

. a!!b.e fraction In equation (31) can be reduced to an ordinary ooiqplex. number 

-and -thd .rccdiv.^ui uf tins noiii'bor is the m^iuilnum a^r^litude of the' steady oscl.lla- 

-tion« The value oon be convertod' to physical units by the formula - - 
■ 

■ Nonitfil 'acceleration per g _ "^^s^ EliLZ— ^ ^ • 
Stick friction in i&ounds PB^o^pg . Z '^^hQQ 

where PhDQ value of elevator dniqping requirod for the condition of neutral 

dynoitdc stability. 
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i'igare 1.- Trim force J-q, ^•nd puH--ap :force as a f-unction of hinge- 
moment parameters. Xa.c. ~ 0.05c: W/Sw = 40 pounds per 
sq-uare foot; c = 7 feet; sea level, i*^ = stick force in pounds for 
.0; 5*21 ^ stick force in pounds per g nonml acceleration* 
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Figure 2.- Trim force ?-a a-ncL pull-up force Jn ^-s a function of hinge- 
moment parametex's, xa,.c. ='0; W/ S'w = 40 poxmds per square 
foot; c - 7 feet; sea level. 2'u = stic^ force in pounds for LY/7 = 
1.0; £"15^ = stick force in pounds per g normal acceleration. 
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Figure 3.^ Trirn force I-c, and pxxll^up force Ftj as a function of hiza^e 

moment parameters. xa,c. - -0,05c; W/s^ ^ 40 pounds per 
square foot; c 7 feet; sea level, ?u = stick force in pounds for 
AV/v - 1.0; ?n ^ stick force in pounds per g normal acceleration. 
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Figure 4.-^ Trim force F-^ ani pull-up force J-^-^ as a faction of hinge- 
moment parameters. Masfe-uji*balanced elevator (h = 10); 
2Ea.c- " O^OSc; W/s^/ = 40 po*arids per square foot; c = 7 feet; sea 
level. Fij. = stici>: force in pouxxds for AY/v = 1.0;? - sticli force 
in pounds per g normal acceleration. 
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Figure 5.^ Effect of altitude on values of stick^force gradient, 
^a.c. ^ 0.05c; w/s^ = 40 pounds per square foot; 

c = 7 feet- 
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Figure 8»- Effect of center-of-gravity location on stat>ility ■boundaries. 
[i, = 13. 5 » 
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9.- JiJffect of elev8tor-control--'System moment of inertia on the 
boundary for increasing oscillations, = 12e5- 
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Jigure 10.- Variation of stability "boLindaries with density parameter,. 
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Figure 11..- Effect of raass unbalance of the elevator control system 
on the stability boundaries, p, = 12.5, 
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Figure 12.- Effect of location of mass unbalance of trie elevator control 
system on the stability 'bo-uncLaries. n = 12.5; ig = 01 ig = 0, 
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Figure IS,- Effect of elevator damping on the bo-undary for increasing oscill- 
ations for various center-of-gravity locations and elevator-mass- 
balance conditions, Cfa^^ = 0.1.; p, - l-£,5» 
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Figure 14,- Effect of elevator damping on the bo"undary for increasing 

oscillations for various center- of -gravity locations and 
elevator-mass-'baiance conditionso Cli^^ " 0*1 j p, = 37*5 
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Figure 15,- Illustration of statilizing effect of elevator damping 
(viscous friction), x^^c. M- 12»5; h = 10. 
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Figure 16*- Illustration of destalDilizing effect of elevator damping 
(viscous friction), ^^^^^^ = -0«05cJ \i = 13.5s = 10* 



•HACA 



Fig. ,17 




-.6 



-.1 



-.5 



-.4 -.3 -.2 

Restoring tendency, 



-.1 



Figure 17.^ Method of constructing boundary for complete damping. 
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Pigure 20.- Tariation of maximim normal acceleration in the steady oscill- 
ation with elevator hinge-moment parameters, x = -.O.OBc* 
^- = 37.5; h - 10; ie = 2; ig = 0. ^-^^ 
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